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Abstract 

We derive a threshold result for fault-tolerant quantum computation for local non-Markovian noise 
models. The role of error amplitude in our analysis is played by the product of the elementary gate time 
to and the spectral width of the interaction Hamiltonian between system and bath. We discuss extensions 
of our model and the applicability of our analysis. 

1 Introduction 

Whether or not quantum computing will become reality will at some point depend on whether we can 
implement quantum computation fault-tolerantly. This would imply that even though the quantum circuitry 
and storage are faulty, it is possible by error-correction to perform errorfree quantum computation for an 
unlimited amount of time while incurring an overhead that is polylogarithmic in time and space, see 0, 
l2l . l4l . 0, 0, (6| and [7]. For this 'software' solution that uses concatenated coding techniques, an 
error probability threshold of the order of 10~ 4 — 10~ 6 per qubit per clock-cycle has been given for the 
simplest error models, meaning that for an error probability below this threshold fault- tolerant quantum 
computation is possible. These estimates heavily depend on error modelling, the efficiency of the error- 
correcting circuits, and the codes that are used. Different and potentially better estimates are possible, see 
for example [ 8 ], [9 1 and [10]. Another solution to the fault-tolerance problem proposed by Kitaev is to make 
the hardware intrinsically fault-tolerant by using topological degrees of freedom such as anyonic excitations 
as qubits I TTI . 

In Refs. and [ 4 ] the threshold result for fault-tolerance is derived for various error models, including 
ones with exponentially decaying correlations. However, this general model of exponentially decaying 
correlations does not make direct contact with a detailed physical model of decoherence. Such a physical 
model of decoherence starts from a Hamiltonian description involving the environmental degrees of freedom 
and the computer 'system' degrees of freedom. 

Starting from such a Hamiltonian picture it was argued in a paper by Alicki et al. [12 ] that fault-tolerant 
quantum computation may not be possible when the environment of the quantum computer has a long-time 
memory. 

In this paper we carry out a detailed threshold analysis for some non-Markovian error models. Our 
findings are not in agreement with the views put forward in the paper by Alicki et al, that is, we can 
derive a threshold result in the non-Markovian regime if we make certain reasonable assumptions about the 
spatial structure and interaction amongst the environments of the qubits. The results of our paper and the 
previous results in the literature are summarized in Section @]of this paper. In section ITTTI we introduce our 
notation and our assumptions on the decoherence model. In section [l"2l we introduce our measure of error or 
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decoherence strength which we motivate with a small example. Then in section [T31 we prove some simple 
lemmas that will be used in the fault-tolerance analysis and in section PHf! we discuss the overall picture of 
a fault-tolerance derivation, in particular the parts of this derivation that do not depend on the decoherence 
model. Then in SectionEJwe fill in the technical details to obtain the threshold result expressed in Theorem 
1. In Section|3]we generalize our decoherence model to incorporate more relaxed conditions on the spatial 
structure of the bath and we discuss further possible extensions. In Section |4] we give an overview of all 
known fault-tolerance results including ours and in the last section|5]we discuss several physical systems in 
which our analysis may be applicable. 



1.1 Notation and Explanation of the Decoherence Model 



We use the following operator norm: \\A\\ = max||^,|| =1 \ \A\ t^}\\ where || 1 = \\ip\\ = \J (ip The 
following properties will be used: \\A + B\\ < \\A\\ + \ \B\\, \\U\\ = 1 if U is unitary, and ||AB|| < 
1 1 .A 1 1 \ \B\\. An operator H that acts on system qubit i or qubits i and j (and potentially another quantum 
system) is denoted as H[qi\ or i?[qj, q 3 -]. A unitary evolution for the time-interval t to t + to is denoted 
as U(t + to, t). to is the time it takes to do an elementary (one or two qubit) gate. The identity operator 
is denoted as I and e denotes the base of the natural logarithm. We will also use the trace-norm denoted 
by \\A\ |i = TrV At A and the classical variation distance between probability distributions P and Q: ||P — 

Qlli = EiTO-Q(OI- 



B 






R 






nS 






A . 





Bl 



ql 



B2 



q2 



(a) 



B1+B2 



ql 
q2 



B3 

q3 






B3 


q3 









t+t 



(b) 

Figure 1: Schematic representation of the model, (a) The system S consists of a register R of qubits plus 
ancillas A that can be reset during the computation. The system S is coupled to the environment, or bath, B. 
(b) The decoherence model. Each qubit q, is coupled to an individual bath Bj. When two qubits interact, 
they may interact with one common bath. 

The following assumptions have been shown to be necessary for fault-tolerance and thus we keep these 
assumptions in our analysis: 

• It is possible to operate gates on different qubits in parallel. 
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• We have fresh ancilla qubits to our disposal. These ancilla qubits are prepared off-line in the exact 
computational state | 00 ... 0) and they can be used in the circuit when necessary. They function as a 
heat-sink which removes entropy from the computation. 

In Figure [^a) three types of quantum systems are sketched that differ in function and in the amount 
of control that we can exert over them. First, there is R, for quantum Registers, that we can control and 
use for our computation. Secondly, there is A, for Ancillas, which are used for error-correction and fault- 
tolerant gate construction during the computation. The systems R and A taken together are denoted as S for 
System of which single qubits are denoted by the letter q. Clean ancilla registers set to | 00 ... 0) are added 
during the computation and can be removed after having interacted with (1) other parts of the system S by 
error-correcting procedures and (2) the bath B according to some fixed interaction Hamiltonian. 

We will assume that the third system, the bath B, which interacts with system and ancillas has a local 
structure, illustrated in Fig.^b). We will generalize this model in SectionPJ] Every qubit (q 1; q 2 , q3-..) of 
the system has its own bath (Bi, B2, B3..). Only during the time when two qubits interact their baths (Bi 
and B2 in the figure) can interact. The idea behind this modelling is that the bath is localized in space, i.e. 
is associated with the place where the qubit is stored. But when qubits interact, they need to be brought 
together and so they may share a common bath. In the picture Bi + B2 at time t + to are. suggested to be the 
same baths that qubits q± and q 2 interacted with at time t, but in general they may also be different baths. 
For example, when qubits and q 2 have to be moved in order to interact, they may see a partially new 
environment at time t + to- This distinction will not be important in our analysis. 

Most importantly, in this model, each bath can have an arbitrarily long memory; at no point in our 
derivation will we make a Markovian assumption. This implies that, for example, the bath Bi may contain 
information about qubit qi at time t, then interact with bath B 2 at time t + to and pass this information on 
to bath B 2 etc. The interaction Hamiltonian of a single qubit of the system (R or A) with the bath is given 
by 

-HsbM = ^Vkfail ® A k- (1) 
k 

with the Pauli-matrices acting on qubit q, and Ap, is some Hermitian operator on the bath of the qubit q, 
which is not equal to the identity I. During a two qubit-gate both qubits may interact with both baths. For 
simplicity (see footnote [ 13 1) we assume that the interaction is of the form 

#sEs[qi, qj] = ^sbM + Hssfaj], (2) 

where the bath part of each H SB [q^] is an operator on the joint bath of qubits q^ and q„-. We do not care about 
the time-evolution of the baths except that it has to obey the "local bath assumption", i.e. noninteracting 
qubits have noninteracting baths. The system (register and ancilla) evolution H RA (t) is time-dependent 
and represents the fault-tolerant quantum circuit that we want to implement. This evolution is built from a 
sequence of one and two qubit gates and, as was said before, to is the time it takes to perform any such gate. 

1.2 Measure of Decoherence Strength 

Our results will depend on the strength of the coupling Hamiltonian ff SB [qj]. There is an additional freedom 
in determining H SB [q^] , namely we can always add a term al s [qj] ® I B where a is an arbitrary real constant 
and I is the identity operator. This is possible since it merely shifts the spectrum (see footnote [ 14)). Let /ij 
be the eigenvalues of H SB . With this freedom we see that 

min||i? SB [qi] + aJ s [q;] <g)I B || = (/i max - Mmin)/2 = A SB [qi], (3) 
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the spectral width of the interaction Hamiltonian (divided by 2). Our analysis will apply to physical systems 
where one can bound 

V qi € S, A S bN < A . (4) 

where Ao is a small constant which will enter the threshold result, Theorem ^ together with to, the funda- 
mental gate time. In what follows we will denote A SB [qj] as A SB or A assuming that the spectral width is 
the same for each qubit in S. 

We justify the use of this norm in the following way. Consider a single qubit coupled to a bath such 
that both bath and system Hamiltonians are zero but there exists nonzero coupling. To what extent will an 
arbitrary initial state of qubit and bath change under this interaction? We can consider the minimum fidelity 
of an initial state V'sb(O) with the evolved state at time t: 

F min (t)=min\(iP(t)\rP(Q))\. (5) 

i/>(0) 

For small times t such that A SB i < 7r/2 the minimum fidelity can be achieved by taking | 0(0)) = 
V'max) + | "0mm)) where I VVax/min) are the eigenvectors of H SB with largest and smallest eigenvalues. 
Then we have 

Fmin(t) = cos(At) « 1 - A 2 t 2 /2 + 0((At) 4 ). (6) 

Note that this fidelity decay includes the effects on the bath. For this reason this fidelity decay overestimates 
the effects of decoherence, in other words F(p s (t), ps(0)) > F m i n . 

One may compare this fidelity decay with that of other decoherence processes, for example the depo- 
larizing channel £ with depolarizing probability p. For such a channel we have F{\ -0)3,5(1 rp)(ip |)s) = 
^jl — | [,.15 1- Thus, loosely speaking, At could be interpreted as an error amplitude whose square is an 
error probability. 

Thus, this brief analysis shows that for some initial states -0 SB (0) the norm of the interaction Hamiltonian 
measures exactly how the state changes due to the interaction. Since our environment is non-Markovian we 
cannot exclude such bad initial states, in other words we cannot assume that the decoherence is just due to 
the interactive evolution of an initially unentangled bath and system. 



1.3 Error Modelling Tools 

The following simple lemma will be used repeatedly in this paper: 

Lemma 1 Let a unitary transformation U = U n . . .U\ where Ui = Gi + Bi and the operator Gi and B{ 
are not necessarily unitary. Let U = B + G where we define B to be the sum of terms which contains at 
least k factors Bi. Let \ \Bi\ \ < e and thus \ \Gi\ \ < 1 + e. We have 

||B|| < ( n k )e k (l + e) n - k . (7) 



If Gi is unitary, we have 



Proof: We can think about U as a binary tree of depth n such that the children of each node are labelled 
with Gi or Bi at depth i. We prune the tree in the following way; when a branch has k factors Bi in its 
path, we terminate this whole branch with the remaining U n . . . U m . The sum of these terminated branches 
is B. B can be bounded by observing that there are (?) terminated branches each of which have norm 
at most | \Bi\ \ k \ \Gi\ \ n ~ k (since each branch is a sequence of Gi transformations interspersed with k Bi 
transformations followed by unitary transformations). □ 
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It is easy to prove the following (see also Ref. Q6)) 

Lemma 2 Consider a time-interval [t, t + to] and a single qubit q £ S which does not interact with any 
other qubit in S at that time. The time-evolution for this qubit is given by some unitary evolution U[q] 
involving its bath B. Let Uo[q] = U s [q] <8> U B be the free uncoupled evolution for this qubit. We can write 

U[q] = U [q] + E[q], (9) 

where E[q] is a fault-operator with norm 

\\E[q}\\ < t \\H SB [q]\\ = t A SB [q] < t X . (10) 



Proof: We drop writing the dependence on qubit q for the proof. For the qubit evolution in the interval, 
using the Trotter expansion we can write 

U= lim U^ =1 (UtU^Ut). (11) 

n^oo 

where U^ 71 is the time-evolution for K = S, B or coupling SB during the time-interval t m of length to/n. 
Now in this expansion we may write Ugg = I — iH SB to/n + O(-fc) and omit these higher order terms. Let 
us call G m = U^Ul™ and B m = -i&Utr H SB Ul m as in Lemma[T] We thus have \ \B m \\ < &\\H SB \\. 
Note that G m is unitary and we have a binary tree of depth n — ► oo and can use Lemma[2with k = 1. This 
gives 

||E|| = ||B|| <i ||#sB||. (12) 

□ 

A similar statement holds when we consider the evolution of two interacting qubits. We have that 

UsB[qi,qj] = U [qi,q.j] +E[qi,qj], (13) 
where ||^[ qi , qi ]|| < 2t A 5B [q] < 2t A . 

1.4 Overall Perspective: good and bad fault-paths 

Since the bath may retain information about the time-evolution and error processes for arbitrary long times 
we cannot describe the decoherence process by sequences of superoperators on the system qubits. Instead, 
there is a single superoperator for the entire computation that is obtained by tracing over the bath at the end 
of the computation. Thus in our analysis we will consider the entire unitary evolution of system, bath and 
ancillas. At time t = bath and ancilla and system are uncoupled and we may always purify the bath, i.e., 
find a pure state in a larger bath Hilbert space which, when the extra Hilbert space is traced out, yields the 
desired mixed state. We can then assume a pure initial product state for the combined system and bath, SB. 
The unitary evolution of the computation consists of a sequence and/or parallel application of the unitary 
gates U[q&, qj](t + to,t) and U[qi](t+to, t). Each such gate, say for two qubits, can be written as a sum of a 
error-free evolution ?7o[<3i) ^'K^ + ^o> t) an d a fault term E[q i: q,-]. Therefore the entire computation can be 
written as a sum over fault-paths, that it, a sum of sequences of unitary error-free operators interspersed with 
fault operators. This is very similar as in the fault-tolerance analysis for Markovian error models, where the 
superoperator during each gate-time to can be expanded in a error-free evolution and an erroneous evolution 
so that the entire superoperator for the circuit is a sum over fault-paths. 
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The main idea behind the threshold result for fault-tolerance is then as follows, see [4|. There are good 
fault-paths with so called sparse numbers of faults which keep being corrected during the computation and 
which lead to (approximately) correct answers of the computation. And there are bad fault-paths which 
contain too many faults to be corrected and imply a crash of the quantum computer. 

Now the goal of our fault-tolerance derivation which is completely analogous in structure as the one in 
141 is to show the following: 

A Sparse fault-paths lead to sparse errors in the computation. This fact uses the formal distinction be- 
tween faults that occur during the computation and the effects of these faults, the errors, that arise 
due the subsequent evolution which can spread the faults. The fact that sparse fault-paths give rise to 
sparse errors is due to fundamental properties of fault-tolerant error-correcting circuitry, namely that 
there exists error-correcting codes and procedures that do not spread faults too much. It is indepen- 
dent of the choice of decoherence model, and can be applied to any model where one can make an 
expansion into fault-paths. See Lemma|3] 

B Sparse errors give good final answers. This is a technical result whose derivation may differ slightly 
in one or the other decoherence model, but which is intuitively sound for all possible decoherence 
models. See Lemma|4] 

C The norm of the operator corresponding to all bad non-sparse fault-paths is "small". This result 
depends crucially on the decoherence model that is chosen, in particular the spatial or temporal cor- 
relations that are allowed. Secondly, it depends on the strength of the errors, that is, only for small 
enough strength below some threshold value will the norm of the bad fault-path operator get small. 
See Section l2~2l 

A,B,C=^ When the bad operator norm is small, the answer of the computation is close to what the good fault- 
path operator yields which is the correct answer according to item B. See Lemma |4] and Theorem 
□ 

Another small comment about our model is the following: In the usual model for error-correction (see 
Ref. [6 1 in FP71 ). measurements are performed to determine the error-syndrome or the correct preparation 
of the ancilla states. Since we prefer to view the entire computation as a unitary process, we may replace 
these measurements by coherent quantum operations. In the error-correction with measurement procedures 
it is assumed that faults can occur in the measurement itself or in the quantum gate that is performed that 
depends on this measurement record, but the measurement record by itself is stable since it is classical. If 
we replace measuring by coherent action for technical reasons in this derivation, it is then fair to assume 
that the qubit that carries the measurement record is errorfree, in other words does no longer interact with a 
bath. This modelling basically allows the standard fault-tolerance results in item A expressed in Lemma|3j 
to carry over in the simplest way to our model. 

2 Threshold result 

2.1 Nomenclature 

Let the basic errorfree quantum circuit denoted by Mq consist of N locations [4|. Each location is given by 
a triple ({q}, G, t) where {q} denotes the qubits (one or two at most) involved in some gate G (G could be 
I) at time t in the quantum circuit. In the following, E[i] or U[i] will denote operators that involve location 
i, i.e. if qi and q 2 interact at location i we will write ?7o[qi, q 2 ] = Uo[i] instead of enumerating the qubits. 
For fault-tolerance one constructs a family of circuits M r by concatenation. That is, we fix a computation 
code C (see definition 15 in Ref. [4|), for example a CSS code, encoding one qubit into (say) m qubits 
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[ 18 1. We obtain the circuit M r by replacing each location in the circuit Mq by a block of encoded qubits to 
which we apply an error-correcting procedure followed by a fault-tolerant implementation of G, see Fig. f2] 
Repeated substitution will gives us a circuit M r at concatenation level r. 



E 




Gfaulttol 




E 





M 



r-1 



Figure 2: Every single or two-qubit gate G in the circuit M r _i gets replaced by an error-correcting procedure 
E followed by a fault-tolerant implementation of G, Gfaulttol (possibly involving ancillas). 



Essential are the following definitions and a lemma taken from Ref. [4] which define sparseness of a set 
of locations and error-spread of a code: 

Definitions from Ref. |4|: 

• A set of qubits in M r is called an s-block if they originate from 1 qubit in M r _ s . A s-working period 
in M r is a time-interval which originates from one time-step in M r _ s . A s-rectangle in M r is a set of 
locations that originate from one location in M r - S . 

• Let B be a set of r -blocks in the circuit M r . An (r, 1) -sparse set of qubits A in B is a set of qubits in 
which for every r-block in B, there is at most one (r — l)-block such that the set A in this block is not 
(r — 1, l)-sparse. A (0, l)-sparse set of qubits in Mq is an empty set of qubits. 

• A set of locations in a r -rectangle is (r, \)-sparse when there is at most 1 (r — l)-rectangle such that 
the set is not (r — 1, l)-sparse in that (r — \)-rectangle. A fault-path in M r is (r, \)-sparse if in each 
r -rectangle, the set of faulty locations is (r, \)-sparse. 

• A computation code C has spread s if one fault which occurs in a particular 1 -rectangle affects at 
most s qubits in each 1 -block at the end of that 1 -rectangle, i.e. causes at most s errors in each 
1 -block. 

• Let Ac be the number of locations in a 1-rectangle for a given code C. 

We state the basic lemma about properties of computation codes which was proved in Ref. El (with a 
correction). 

Lemma 3 (A: Lemma 8 in |4| with a correction) Let C be a computation code that can correct 2 errors 
and has spread s = 1. Consider a computation M r subjected to a (r, 1) -sparse fault-path. At the end of 
each r -working period the set of errors is (r, 1) -sparse. 

Thus for simplicity we will be using a quantum computation code that encodes one qubit and can correct 
two errors and has spread s = 1. We denote the entire unitary evolution of M r including the bath as Q r . We 
may write Q r = Q r G + Q r B where Q T G is a sum over good (r, 1) -sparse fault-path operators and Q T B contains 
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the bad non-sparse terms. A fault-path operator E SB that is (r, l)-sparse is a sequence of free evolutions 
Uo[i) for all locations except that in every r-rectangle there is a (r, l)-sparse set of locations where a fault 
operator E[i] occurs. 

Definition 1 (Operators in the Interaction Picture) Let Uo(t 2 ,ti) = U s (t2,ti) (g> U B (t2, h) be the free 
uncoupled evolution of system and bath in the time-interval [ii, £2]- We define a fault-operator E(t2,ti) in 
the interaction picture as 

E(t2,h) = UofatjEUlfah). (14) 

The interpretation is that E(t2,t\) is the spread of a fault E that occurs at t\ due to the subsequent free 
evolution. 

Then it is simple to see the following: 

Proposition 1 (Error Spread in the Interaction Picture) Consider a quantum circuit M. Let Uo(tp,tj) 
be the uncoupled evolution for M. Faults occur at a set of 'time-resolved' locations 

T = ((h,h), («2,*2), • • • , (ik,tk)) 

where i\, . . . , if. is the set of distinct locations of the faults and t\, . . . , ifc label the specific times that the 
faults occur at the locations. Let E SB (T) be a particular fault-path operator in which at every faulty 
location (i, t) £ T we replace Uq[i] by a fault-operator E[i]. We have 

EgB^uiitF,^) = E[i k ](t F , t k ) . . . E[h](t F , h). (15) 

We note that the system part of E sb Uq is I everywhere except for the qubits that are in the causal cone of 
the faulty locations, i.e. the qubits to which the errors potentially have spread. 

Proof: This can be shown by inserting I = U^tp ,ti)Uo(tp ,ti) in the appropriate places and then using 
the definition of fault operators in the interaction picture. □ 

Now we include error-correction and differentiate between the ancilla systems A used for error-correction 
which may contain noise and the registers R in which the errors remain sparse. Note that all these ancillas 
are in principle discarded after being used, but we may as well leave them around. Let K\c be the restriction 
of the operator K to vectors in the code-space of C, i.e. K\c = KPc where Pc* is the projector on the 
codespace. 

Let us consider a fault-path operator E SB representing a single fault E at time t on some block that 
is subsequently corrected by an errorfree error-correcting procedure. Let | IN) be the initial state of the 
computer, bath and ancillas and Uo(tF,ti) be the perfect evolution. We have 

E SB \ IN) = E SB U j U \ IN) = E SB ut\ MtF))- (16) 

where | V'c(^f)) is the final perfect state of the computer, prior to decoding and therefore in the code-space. 
E SB is the sequence Uoftp, t)E Uo(t, to) where Uo(tF,t) includes the error correction operation. In other 
words, in the interaction picture, we can write 

E SB \W) =E(t F ,t)\i/> c (t F )), (17) 

The error-correcting conditions (see [ 15 1, par. 10.3) imply that when acting on the code space and an ancilla 
state set to | 00 ... 0) the operator E(tp, t) will be E(t F , t) = l\c (Junk) AB where Junk is some arbitrary 
operator on the ancilla (that receives the error syndrome in the error-correcting procedure) and bath. In Eq. 



2 THRESHOLD RES ULT 



9 



(flTt the final errorfree state has all ancillas set to | 00 ... 0) and the system state is in the code-space and 
thus the error acts as I on the system. 

Similarly, let E SB contain two faults at times t% < t 2 that have not spread (say) and are then corrected 
by a perfect error-correcting procedure. We have 



Let us assume, for example, that E\ occurs prior to error-correction and E 2 occurs during error-correction. 
Then due to the error correction Ei(tF,ti) acts as I on the code space when the ancilla used for error- 
correction is set to I 00 ... 0) and acts as Junk on this ancilla and the bath. The error E2 will not be 
corrected and may still be present (but will not have spread to more qubits in the block due to the spread 
properties of the code that is used) after error- correction. Thus in total we can write for this process that 
Ei(tf,ti) acts as I on the code space, whereas £2(^7 ^2) is an operator that acts on the code space as at 
most one error per block. 

Alternatively, both faults could occur prior to error-correcting so they can both be corrected by our code. 
This implies that both Ei(tF,h) and £2(^7 ^2) act as I on the code-space. Note that after the first fault the 
ancilla will be partially filled (i.e. not be | 00 . . . 0)) but since the code can correct two errors there is still 
space to put the second error syndrome in. However a third operator ^(i^, £3) would no longer act as I on 
the code-space since the code cannot correct three errors. 

In other words, with these examples we can see how Lemma |3] can be translated in terms of the sparse- 
ness of the errors in the interaction picture, i.e. the sparseness of places where they act as non-identity 
on the final encoded state of the register qubits. In the next lemma we need to consider the effect of such 
sparse fault-path operators E SB on the final state of the computer. This is the state of the computer obtained 
after fault-tolerant decoding which is as follows. The fault-tolerant decoding procedure for a single level of 
encoding takes a codeword | c) and 'copies' (by doing CNOT gates) the codeword m times. Then on each 
'copy' we determine what state it encodes and then we take the majority of the m answers. This procedure 
is done recursively when more levels of encoding are used. 

In the fault-tolerant decoding procedure faults can occur on the codewords, i.e. as incoming faults, 
during the copying procedures and during the determination of what is encoded by the codeword. The last 
procedure will usually be a conversion from a quantum state to a classical bit string since this will be the 
most efficient. This implies that the step of taking the recursive majority of these bits is basically errorfree 
since it only involves classical data. In the next Lemma we model this by coherent quantum operations that 
output superpositions of decoded bit strings followed by an error-free measurement that takes the recursive 
majority of these bits. 

Lemma 4 (B: Sparse faults give almost correct answers) Let Q r = Q r G + Q r B the unitary evolution of 
M r and let \ \Q r B \\ < e < 1/2. Let Po(i) be the output probability distribution under measurement of some 
set of qubits of the error-free original computation Mq. Let P(i) be the simulated output distribution of the 
encoded computation M r with evolution Q r . We have 



Proof: The initial state of the computer is | IN) RAB = | 00 . . . 0)ra <8> | IN B ) for some arbitrary state | IN B ). 
Let Uq be the error-free evolution of M r including the final decoding operation. Thus let Uq\ IN) RAB = 
I OUT ) R ® I REST) AB . Let Q r \ IN) RAB = | OUT) RAB and Q r B/G \ IN) RAB = | OUT B/G ) RAB . We will drop 
the label RAB from now on. The norm of | OUTc) will be denoted as | |OUTg| |. We have 



E SB \ IN) = E 2 {t F , t 2 )E x (t F ,t x )\ 1>c(t F )). 



(18) 



||P -P||i < \/2l + 16e. 



(19) 



l = ||Q r |IN)||<||%|IN)|| + ||Q!B|IN)|| 



(20) 
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so that ||OUT G || > 1 - ||<%|IN)|| > 1 - e. On the other hand \\Q r G \\ = \ \Q r - Q r B \ \ < 1 + e. 

Let G be the set of (r, l)-sparse fault-paths. We have Q r G = ^2q- e c E SB (T) where E SB (T) is the 
fault-path operator of a (r, 1) -sparse fault-path labelled by location and time index set T. We can write 

| OUT G ) = Y, E SB (T)U$\ OUT ) R ® | REST) AB . (21) 

TeG 

By the arguments above and the fundamental Lemma[3]we know that E SB (T)U^ is I everywhere except on 
a (r, 1) -sparse set of qubits. Let w be the number of output qubits of Mq. The ideal state | OUTo)r has the 
property that all qubits in an r-block have the same value in the computational basis, i.e. 

|OUT ) R = £ a n .., w \hr mr ...\i w r m \ (22) 

where m is the number of qubits in a 1-block. The final step of the computation is a measurement of 
all output qubits that takes the recursive majority on the block to get the final output string i of length w 
with probability ¥ tot (i). We model this measurement using POVM elements E k , - J2k-^k = I- Since 
not all these w output bits may be relevant output bits of Mq, we may use the fact that trace-distance is 
non-increasing over tracing lU31 so that 

||Po < \K ot -lP to *lli, (23) 

where F tot (k) = TrE k \ OUT)(OUT (rab and P£°*(ifc) = TiE k \ OUT )(OUT | R . Let us also define 
Pg*, the distribution of outcomes if the state of the computer would be the normalized state | OUTq) = 
| OUTg)/||OUTg||. The triangle inequality and the properties of the trace-norm imply that 

| |P*°* - IP* oi 1 1 x < ||P to * - P(j*||l + ||Pg* - IPo 0< 1 1 1 < 
|| | OUT) (OUT | - | OUTg)(OUTg | ||i + ||P(f - P o *||i- (24) 

Here the first term can be bounded, using the relation of the trace norm to the fidelity F(ip, 4>) = \{ip\(f)}\ 
O, as 



OUT) (OUT | - | OUT^) (OUT^ | 1 1 1 < J I - F (OUT, OUT^) 2 < 



y/l - ||OUT G || 2 < y/2e - e 2 . (25) 

Now consider the second tracenorm on the r.h.s. of Eq. <!24b . We note that all states that are linear combina- 
tions of (r, 1) -sparse error sets applied to the state | ki)® mr . . . | k w )® mT will give rise to the measurement 
outcome k since we are taking majorities. We can model E k = P k where P k is the projector onto the space 
of computational basis states that give rise to the majority output string k. Thus we have 

P k \ OUT G ) = a kl ... kw Esb(T)U^\ hf mr ... | k w f mr | REST) AB . (26) 
TeG 

which can be written as oc kl .., kw Q G \ ^k) f° r some normalized state | ifik) rab- This implies that the second 
term in Eq. d24l can be bounded as 



^-p °*||i = 5> fc ,..U 2 



IIQ G I^)II 2 1 

IIOUTolP ' 



< 



,2 tj 

k 



^ \a kl ...k w \ max 



\\Qh\^)\\ 2 

||OUT G ||2 



* (T^P (27) 



using the bounding inequalities of HQgII an d ||OUT G ||. All bounds put together, using e < 1/2, give the 
result, Eq. ©. □ 
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2.2 Step C: non-sparse fault-paths have small norm 

Consider the evolution Q r which can be viewed as a sequence of unitary evolutions, one for each r- 
rectangle, since qubits in different rectangles do not interact. The number of locations in Mq is N. The 
computation Q r is bad when at least one r-rectangle is bad, or using LemmaU 

\\Q r B \\<N\\R r B \\\\K" G \\ N -\ (28) 

where R r B and R! G are the good and bad parts of the unitary evolution R r for some r-rectangle. The unitarity 
of R r implies that we can bound 1 1 R r G \\ < 1 + 1 1 R r B 1 1 . In each rectangle we can view the entire evolution 
as a sequence of unitary evolutions, one for each (r — l)-rectangle. Note that we are again using the fact 
that non-interacting qubits have non-interacting baths. A r-rectangle is bad when there are at least two 
(r — l)-rectangles which contain sets of faulty locations which are not (r — 1, 1) sparse. This implies, using 
Lemma [J again, that 



II^II^^Jll^- 1 !! 2 !!^- 1 !!^- 2 , (29) 

where we can use that 1 1 R r G 1 1 1 < 1 + 1 1 R r B ~ 1 1 1 . When r = 1, R G is a unitary operation and thus 1 1 R G \ \ = 1. 
This recurrence in r is identical to the one in Lemma 1 1 in Ref. [4 1 and thus the solution and results are the 
same if we replace rj in Ref. |4j by Aoio- Thus the critical error threshold value is 

^ = eAcdc-ty (30) 

Here we can observe a difference with the simplest error model with error probability p for which the 

critical value is p c = ,} . [4|. The dimensionless quantity Aoio plays the role of an amplitude, see Sec. [5] 

( c i 



2 



which implies that this threshold value may be more stringent than in the simple probabilistic error model 
(see also the critique by Alicki [ 19 1 on our results). However, we believe that this analysis is too course to 
really give information about the value of the threshold. The fact is that in practice, baths do not have infinite 
memory times since they are coupled to many other degrees of freedom. Representing the coupling between 
bath and system as a pure coherent evolution was needed in this analysis to deal with the non-Markovian 
dynamics; however, we do not expect this formal procedure to give rise to an optimal error threshold. 

The idea of the remaining derivation given in Ref. is to show that when Ao^o < (Aoio)e for large 
enough concatenation level r, | \Q r B \ | < e for arbitrary small e. Lemma|4]then tells how much our quantum 
computation errs from the error-free computation. Summarizing we get the following, as in Ref. |U: 

Theorem 1 (Threshold Theorem for Local Non-Markovian Noise) Let N be the number of locations of 
an errorfree quantum computation M that outputs samples from a probability distribution P. There exists 
a quantum circuit M' subjected to noise according to the Hamiltonian H SB and bath Hamiltonian H B 
that obeys the "noninteracting qubits have noninteracting baths" assumption which outputs the probability 
distribution P' such that 

||P'-P||l<e, (31) 
when A SB to < e A c (A c ~i) an< ^ ^ c iJ ^ e num ^ er of locations in a 1-rectangle of a (computation) code C 
that can correct two errors and has spread s = 1. The number of locations in M' is A r polylog-^. 

Proof: This follows from the basic results in Ref. @ and Lemma |4] and step C above. It was proved in 
Ref. [4] for Ao^o below the threshold (Ao^o)c hi Eq. d30l when the concatenation level r = c\ log(log y + 
C2) + C3 for constants c\, C2 and C3 we have HQ^H < e'. So we choose M' = M r , the computation at this 
concatenation level r which implies that ||P' — P||i < + 16e' = e. The number of resources (time 
and space, related to Ac) in M' scales exponentially, i.e. the number of locations in M r is NA G . With the 
dependence of r on N and e' this implies the polylogarithmic overhead in terms of N and yfe. 

□ 
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3 Extension to Decoherence Models with Clustered Qubits 

In the most general noise model we start with a Hamiltonian description of system and bath. We will assume 
that such Hamiltonians are l-systemlocal, that is, the interaction Hamiltonian between system and bath is a 
sum of terms each of which couples a single qubit to some part of the bath. This covers many interaction 
Hamiltonians in systems that are being considered for quantum computation (see footnote l20l l. 

We have seen that basically the only place where the noise model enters the derivation of fault-tolerance 
is in section l2~2l i.e. the derivation that the total amplitude/probability/norm for non-sparse fault-paths at 
concatenation level r goes (doubly exponentially fast in r) to zero when the initial error strength is below 
the threshold. Locality of the interaction Hamiltonian is an important (and necessary) ingredient in the 
derivation of fault-tolerance since it implies -without any further assumptions on the structure of the bath or 
the (Non)-Markovian character of the system- that fault-path operators with k faults have a norm bounded 
by (2Ao£o) fc ( see Appendix [Al. This bound is not strong enough by itself to derive that HQs 1 1 becomes 
arbitrarily small for sufficiently small Aoio- We find that there are technical and potentially fundamental 
problems in the derivation of step C, for the most general local Hamiltonian model both in the Markovian 
case as well as in the non-Markovian case. The problems are due to the fact that all qubits of the computer 
potentially couple at a given time to the same bath which was prevented in the derivation of Theorem 1 by 
assuming that "noninteracting qubits have noninteracting baths". The problem is basically due to the fact 
that the unitary evolution of a working period cannot be written as a product of unitary evolutions for each 
rectangle in the working period since different rectangles may share their bath. 

We thus need to consider restricted models that are still physically very relevant: 

3.1 Clustered Qubits at encoding level r = 1 

We can generalize the model in Sec. 11.11 i.e. noninteracting qubits have noninteracting baths, to one in 
which a cluster of qubits can share a bath. The model is depicted in Fig. |3] We will assume that qubits that 
are contained in a 1 -rectangle of M r may share a common bath whereas qubits in different 1 -rectangles do 
not share a bath. We imagine that baths are attached to physical locations, so that the interaction regions of 
different r-rectangles are physically separate. This means that from one 1-working period to the next one, 
qubits have to be moved around, i.e. qubits that participate in one 1-rectangle have to be brought together. 

Let us for the moment neglect the machinery that is necessary to move qubits around. Then we can 
observe that the entire computation M r can be viewed as a sequence of unitary gates each involving a single 
1-rectangle. In the 1-rectangle we cannot decompose the evolution as a sequence of unitary transformations 
for each location. But at this lowest level r = 1 it is simple to derive a bound on the bad part R l B of the 
unitary operation R l . Given this bound we can insert it in the previous recurrence of Eq. d29t and determine 
a threshold which is the same as before. Here is the bound on R B : 

Lemma 5 Let R 1 = R l B + Rq be the unitary transformation of a 1-rectangle where R B is a sum of 
non-sparse fault-path operators, i.e. each such operator contains at least two locations with faults. Then 

\\R B \\ <2(A c t \ ) 2 . (32) 

and \\RqW < 1 + \\R B \\- 

Proof: We do a Trotter expansion for R 1 as in Lemma[6]and obtain a tree with infinite depth. We combine 
branches of the tree in the following way: 1. after a location has become faulty we append the full unitary 
for the remaining time of the location and 2. if two faults have occurred at two different locations we do 
no longer branch the tree and just append the entire remaining unitary transformation to that branch. In this 
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Bath of cluster 1 




Bath of cluster 1+2 



3ath of cluster 2 



Gfaulttol2: 




[t,t+A ] 



[t,t+2t x ] 



Figure 3: Schematic representation of decoherence model where clusters of qubits can share a common 
bath. Logical qubits 1 and 2 are encoded once in a block of qubits. In the original circuit these qubits 
first undergo single qubit gates G\ and G 2 and then interact in G3. In encoded form this implies three 
1 -rectangles that each take some time t±; these are denoted by the boxes with fat lines in the figure. Each 
1 -rectangle or cluster has its own bath. These baths may change over time, that is, the bath of cluster 1 may 
evolve or change and not be the same as the environment that this block of qubits sees later. 



way the norm of every time-resolved branch with at least two faults E 2 + (T) is bounded by | \E 2 + (T)\\ < 
\2 tw„ „ ra (Acton/t^ such baches an d thus 



(2Aoio/fi) 2 - There are 



I^bII < (2A ( 



t/nf( Act ° n/t 



< 2{A c t X ) 



(33) 
□ 



A physical example of this decoherence model is the proposal for scalable ion-trap computation 1231 . 
A few qubits are stored in an ion-trap where they may share a common bath. The states of qubits can be 
moved around to let them interact. A small cluster of ion-traps may be used to carry out the fault-tolerant 
circuits and error-correcting at level r = 1 of encoding. 

The issue of moving qubit states around is not entirely trivial and will be addressed in detail in a future 
paper l24l . 



4 Overview 



We would like to summarize the known results, including the ones in this paper, on threshold results for 
different decoherence models. The simplest model is one in which we assume that each location undergoes 
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Spatial Correlations 







Single Location 
Baths 


Cluster Location Baths 
at r = 1 


Arbitrary 
Baths 


Temporal 
Correlations 


Markovian within 
gate-time to 


/ 


/ 


? 


Non-Markovian with 
finite memory time t > to 


/ 


/ 


? 


Non-Markovian with 
unlimited memory time 


/ 


/ 


? 



Table 1 : A / indicates that a fault-tolerance result exists whereas a question mark ? indicates that it is 
not known to exist so far (neither has it been disproved). The results for non-Markovian baths assume 
a 1-systemlocal interaction Hamiltonian that can be bounded in norm. They also assume that we can do 
two-qubit gates between any two qubits in the circuit (that is, we do not take physical locality constraints 
into account). The assumptions on the structure of the system-bath interaction and the bath Hamiltonian are 
given by the three columns. Single Location Baths implies that the interaction and the baths are constrained 
so that for each elementary time-interval (clock-cycle) [t, t + to] the following condition is obeyed: qubits 
that do not interact can only interact with baths which do not interact, see Fig. ^b)- Note that the particular 
baths with which the qubits interact may change over time. Cluster Location Baths is the extension of this 
model covered in Section|3]where a cluster of qubits can share the same bath, see Fig. |3] In the last column 
there is no constraint on the bath. 

an eiTor with probability p and undergoes no error with probability 1 — p. This is a specific example of a 
Markovian model in which in every location has its own separate environment, i.e we have "Single Location 
Baths", see the upper left entry in the Table 1. Generalizations of this model exist 61 0; in these models a 
superoperator S(p) = So(p) +£{p) where Sq corresponds to the error-free evolution and £ to the erroneous 
part, is associated with each location. Again this corresponds to the upper left entry in the table. This model 
has been generalized to allow for more general correlations in space and time in the following manner. In 
Ref. [4] fault-tolerance was derived in a model where it is assumed that the probability for a fault-path 
with k faults is bounded by Cp k (l — p) N ~ k where N is the total number of locations in the circuit (note 
the difference with Eq. (I43t in Appendix lA"b. Similarly, in Ref. [3| fault-tolerance was derived under the 
assumption that a fault-path with at least k faults has probability bounded by Cp k for some constant C. Let 
us call these conditions the exponential decay conditions. Note in the Table that it is not known whether one 
can derive fault-tolerance for a entirely Markovian model but with extended spatial correlations between the 
baths, i.e. for every clock-cycle we have a superoperator that acts on all qubits of the system; the point is 
that it is not clear whether such a superoperator would obey some sort of exponential decay conditions. 

5 Measures of Coupling Strength and Decoherence 

In our analysis the role of error amplitude is played by the dimensionless number Aoto which captures the 
relative strength of the interaction Hamiltonian as compared to the system Hamiltonian. It is this quantity, 
Aoto> that should be O(10~ 4 ) as was determined for some codes. In a purely Markovian analysis we 
typically replace Ao by an inverse T2 or T\ time and this may give a more optimistic idea of the regime of 
fault-tolerance. Let us consider a few examples of decoherence mechanisms and see how sensible it is to use 
A sb as a bound for decoherence. A good example of a non-Markovian decoherence mechanism is a small 
finite dimensional environment localized in space, for example a set of spins nearby the system of interest. 
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An example is the decoherence in NMR due to interactions with nuclear spins in the same molecule. In 
NMR the nuclear exchange coupling between spins a and b is given 

Hint = Jab la ■ h- (34) 

If the J-coupling is treated as a source of decoherence as compared to the Zeeman-splitting ujq for an 
individual spin, then J /ojq can be ~ 1(T 6 (see footnote 11271 1. 

For some physical systems a source of decoherence is a bath of spins, each of which couples to a 
single qubit. An example is the electron spin qubit in a single quantum dot which couples via the hyperfine 
coupling to a large set of nuclear spins in the semiconductor |2"81 . The interaction Hamiltonian is as follows 

TV 

Hint = ^2 a ^ ° ' ^ 35 ) 

i=l 

where a(i) = Avo\tp s (i)\ 2 and A is the hyperfine coupling constant, v the volume of the crystal cell, and 
^(i)! 2 is the probability of the electron to be at the position of nuclear spin i. If we bound J2i | "0 S (^) 1 2 < 1 
we have that ||ff int || < CAvq where C is a small constant (of order 1). This may give a somewhat weak 
upper bound on the decoherence, since we are basically adding the effects of each nuclear spin separately. 

A third type of decoherence mechanism exists which is essentially troublesome in our analysis. This is 
the example of a single qubit, or spin, coupled to a bosonic bath. The interaction Hamiltonian is that of the 
spin-boson model [ 29 1 

JV 

Hint = <Tz® y^X^i + c i a i)> ( 36 ) 

2=1 

where i labels the ith bosonic mode characterized by frequency Wj. The ith bosonic mode has Hamiltonian 
H Ei = u>i(a\ai + 5). In order to represent a continuous bath spectrum, one lets N go to infinity. In that 
limit the coupling constants Cj are determined by the spectral density J(uS) = \ci\ 2 S(u> — ui). The 
spectral density can have various forms, matching the phenomenology of the particular physical system, 
an example is the Ohmic form in which J(co) = a.uie~ U} / ulc where a is a weak coupling constant that has 
physical relevance and lo c is a cutoff frequency that is also determined by the physics. It is clear that 1 1 Hi n t \ \ 
has no physical meaning since it is infinite, the reason being that there are infinitely excited bath states with 
infinitely high energy. We can determine an energy-dependent upper bound on this norm; using properties 
of the norm, we can estimate 

1 1 -Hint I ^>sb|| = 1 1 ^ CiCii + c*a\)\ V>)sb|| < ^2 l Ci l(ll ai l ^)sb|| + ||4l ^)sb||) 

i i 

<J2^T= V^I-HbJ^sb- (37) 



Using the Schwartz inequality we get 



ffintlVOsBll <2 J^^(^\H B \^) SB = 2J (H B )^ SB J™ du^. (38) 



for some state of system and bath | ^} SB where the bath Hamiltonian H B = H Bi . The idea is that for the 
physically relevant states of the bath (H B )^ is bounded. The problem remains that this bound will in general 
be too poor to be physically relevant, since this energy bound may be quite large. Also, for ohmic coupling 
(for example) we have the integral J °° dw^^- = a u c , i.e. linear in ui c . The cutoff uj c may be quite large 
and it is more typical to see decoherence rates depend on log lo c as in the non-Markovian analysis of Ref. 
[ 30 1 for example. 
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5.1 Cooling assumption 

Some progress can be made in finding good bounds for in the case of a bosonic environment if 

additional assumptions about its state can be made. What is troublesome about the potential nonequilibrium 
state of the bath is that expectation values such as TrajOjl ip}(tp |sb and Traf\ ip)(ip | S b mav not be zero 
since the bath state may not be diagonal in the energy or boson number basis. On the other hand, interaction 
with other environments, for example by means of cooling, will dephase the state of the bath (due to energy 
exchange) and drive it to a state that is diagonal in the energy eigenbasis. Under that assumption only the 
terms Tr a\a.- t \ ip)(ip |sb and Tr a;aj| tp)(tp |sb are nonzero. In that scenario, Eq. d37l simplifies to 



H^I^sbII ~^|q| 2 (^ |ata; + 1/2)| V)sb- (39) 

Still, H-ffjntl -0) sb 1 1 can be very large if some modes of the environment are highly excited, rn = a|a; 3> 1. 
However, in a realistic setting, this will be prevented by cooling the bath, i.e. by constantly removing energy 
from it. Without making a Markov approximation, we can thus assume that the occupation numbers n; of 
the bath are upper-bounded by those of a thermal distribution with an effective maximal temperature T e g. 
This gives 

\\H int \iP) SB \\ < JJ^ J(u;)coth(/? cfr cj/2)/2, (40) 
where f3 e g = k J Ts ■ We can evaluate Eq. d40T) for the Ohmic case using Mathematica 




l^n^> SB ||;< A /^/-u; c 2 + — (41) 



where ^'(x) is the first derivative of the digamma function fy(x) = where T(x) is the Gamma func- 



ft 
IX7T 



tion. For „ ^ <C 1 we do a series expansion and obtain 

l|^int|^)sB||^4/^W^ + -^-(^- + Ol-^) ). (42) 




Unlike Eq. d38b . this bound does not involve extensive quantities, such as the total energy (H B )^ SB of the 
bath. However, Eq. (l4*2l still involves the high-frequency cut-off u c because of the zero-point fluctuations 
of the bath. 



6 Conclusion 

Some important open questions remain in the area of fault-tolerant quantum computation. Most importantly, 
is there a threshold result for non-Markovian error models with system-local Hamiltonians and no further 
assumptions on the bath? Is this a technical problem, i.e. how can one efficiently estimate Q' B , or are 
there specific malicious system-bath Hamiltonians that have such effect that the norm of the bad faults does 
not become smaller when increasing r? The next question is whether a better analysis is possible for the 
spin-boson model, which is a highly relevant decoherence model. One would like to evaluate the effect of 
H SB in the sector of physical states, but the characterization of these physical states is unclear due to the 
non-Markovian dynamics. For real systems one is probably interested in a finite memory time r > to which 
may be simpler to solve. For example, one can derive the superoperator for a single spin qubit coupled to 
a bosonic bath in the Born approximation [ 30 1, however a derivation involving more than one system qubit 
may be too hard to do analytically. 
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A Bounds on Fault-Path Norms 

Lemma 6 (Fault-Path Norms) Consider the entire unitary evolution Q r of a quantum computation on SB. 
Let Vq G S, A SB [q] < Xq. We expand Q r as a sum over fault-paths which are characterized by a set of 
faulty locations I. A fault-path operator with k faults has norm bounded by 

\\E{l k )\\<(2\ t ) k . (43) 



Proof: We do a Trotter-expansion of Q r and obtain a tree with an infinite number of branches each of which 
corresponds to a certain time-resolved fault-path. Every time a fault occurs at some time t m and location i m 
we append unitary evolutions for the remaining time of the location, since we do not care that more faults 
occur in that time-interval, the location has failed anyway. These time-resolved fault-path are characterized 
by an index set T = ((ij, t\), (^2,^2), • • • > where ii, . . . , i k is the set of locations of the faults 

and t\, . . . , tk label the specific times that the faults occur at the locations. Every such time-resolved fault- 
operator with k faults has norm bounded 

\\E(T k )\\ ^(^r)"- ( 44 ) 

Now we need to group these time-resolved fault-paths corresponding to faults at sets of locations. For fixed 
n faults can occur in time-intervals of length t/n and thus during a time to ^ time-resolved faults can 
occur. This implies that 

\\E(I k )\\ < || Yl \\E(T k )\\<(2X t ) k . (45) 

□ 
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